Interplay of Dirac fermions and heavy quasiparticles in solids 
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Many-body interactions in crystalline solids can be conveniently described in terms of quasiparticles with 
strongly renormalized masses as compared to those of non-interacting particles. Examples of extreme mass 
renormalization are on the one hand graphene, where the charge carriers obey the linear dispersion relation 
of massless Dirac fermions, and on the other hand heavy-fermion materials where the effective electron mass 
approaches the mass of a proton. Here we show that both extremes, Dirac fermions like they are found in 
graphene and extremely-heavy quasiparticles characteristic for Kondo materials, may not only coexist in a solid 
but can undergo strong mutual interactions. Using the example of EuRli2Si2 we explicitly demonstrate that 
these interactions can take place at the surface and in the bulk. The presence of the linear dispersion is imposed 
solely by the crystal symmetry while the existence of heavy quasiparticles is caused by the localized nature of 
the 4f states. 

PACS numbers: 79.60.-i, 71.27.+a, 74.25.Jb 



The idea that collective excitations, e.g. lattice vibrations, 
and many-body phenomena, e.g. electron correlation, can be 
described in terms of a rescaled single-particle problem d 
was a major breakthrough in condensed matter physics. The 
solutions of such models are quasiparticles which relate the 
many-body interactions to the properties of a virtual single 
particle. This concept can provide an understanding of com- 
plex phenomena like e. g. the exotic low-temperature physics 
of heavy-fermion systemsintermetallics synthesized on the 
basis of 4/ or 5/ elements where the conduction electrons 
act as particles with huge masses, comparable to that of a 
proton OO. The large effective masses are reflected in the 
low temperature properties of the material, for instance, by 
large values of the specific heat or an increase of electric re- 
sistivity below the so-called Kondo temperature. In angle- 
resolved photoemission (ARPES) data, the respective states 
can be identified as flat running states with an almost com- 
plete absence of dispersion. A diametrically-opposed exam- 
ple are charge carriers in graphene, that behave as massless 
Dirac fermions (3H8), which are, like photons, characterized 
by a linear dispersion relation [ 6 ] . The results of A. Geim and 
K. Novoselov explicitly demonstrate a correlation of this lin- 
ear dispersion with other unique properties of graphene like 
the large opacity, the huge charge carrier mobility, and the 
anomalous quantum-Hall effect (U [71 [8). Here we report 
that both extremes, massless and ultra-heavy fermionic quasi- 
particles may not only coexist in a solid but even exhibit mu- 
tual interaction and form composite quasiparticles. 



To this end, we studied the electronic structure of the 
layered material EuRh 2 Si2 by means of ARPES and den- 
sity functional theory (DFT) calculations. EuRh 2 Si2 is a 
stable-valent representative of the broad family of famous 
heavy-fermion compounds like YbRh 2 Si2, CeCu2Si2 and 
URu2Si2 0HO). Around the T-point of the Brillouin zone 
(BZ) we observe a sharp, linearly dispersing band of Rh 4d- 
derived states. Its conical shape is directed towards the Fermi 
level (Ep) and it passes through an Eu 4f 6 ( T Fj) excited 
state multiplet with whom it interacts. A similar interaction 
is maintained when going from the surface into the bulk as 
proven by using soft X-ray ARPES with high bulk sensitiv- 
ity. The discovery of quasiparticles simultaneously formed 
by massless and heavy-fermion states is quite remarkable be- 
cause species showing either Dirac or heavy fermions typi- 
cally exhibit opposing material properties. 

The linear band dispersion around T seems to exist in 
various members of the tetragonal body-centered ThCr2Si2 
(122) family. In the isostructural Kondo lattice compound 
YbRh 2 Si2, a respective band has been previously found but 
there the apex of the Dirac-like cone is slightly deformed by 
hybridization with 4f states at the Fermi level |[TQllT2ll . A 
similar interplay between heavy 4f and massless 4d quasipar- 
ticles also occurs in EuRri2Si2 leading to a pronounced non- 
crossing behaviour that is explicitly reminiscent to that of po- 
lariton generation in optics where a light cone interacts with 
an almost non-dispersive optical phonon branch fT3ll . In con- 
trast to YbRh 2 Si2 the linear band dispersion is preserved at 
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FIG. 1 : Experimental ARPES spectra taken along the T-X-direction 
with the photon energy 120 eV, i.e. sensitive to the Eu 4f emission 
while the Rh 4d states are strongly suppressed by the Cooper min- 
imum, for (a) Si terminated and (b) Eu terminated surfaces. The 
simulated model ARPES spectra for the same surfaces are presented 
in (c) and (d), respectively. 



Ey in EuRh2Si2 due to the absence of energetically low-lying 
4f excitations. Because of its hybridization with the 4/ multi- 
plet the linearly dispersing band accumulates a certain weight 
of 4f character which it carries all the way up to the Fermi 
level region. A computational analysis provides a hint to spin 
polarization of the linear band. This is in close analogy to 
what is observed for topological insulators fT4ifT6l where re- 
spective linear dispersive surface states bridge the band gap 
and lead to metallic behaviour in the surface layer. How- 
ever, the linear state in EuRh 2 Si2 is probably not protected 
by time-reversal symmetry like in topological insulators since 
EuRli2Si2 is a metal and the state is also found in the bulk. 
Thus it may be closer to the case of graphene where the linear 
dispersion emerges due to the lattice symmetry. 

EuRli2Si2 usually cleaves along the Eu layers leaving be- 
hind surfaces terminated either by Si or Eu atoms. In figs.[lj 
and b we compare the electron band structure as seen by 
ARPES along the T-X-direction for the two possible surface 
terminations after cleaving. The band topologies look rather 
different in both cases: The 4f electron emission from bulk 
Eu atoms, visible as the narrow, string-like 4f 6 final state 
multiplet between 0.2 and 0.8 eV binding energies (BE), is 
detected for both surfaces, while the presence of Eu atoms at 
the surface gives rise to an additional broad 4f signal between 
1.2 to 2 eV BE. 

The strong localization of 4f electrons leads to the ob- 
served almost dispersion-less behaviour corresponding to the 
limit of free ions or - rephrased in the band structure picture 
- infinitely heavy quasiparticles. These "horizontal bands are 
disturbed at those points in &-space, where they are crossed 
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FIG. 2: LDA based band structure calculations for a surface covered 
by (a) Si and (b) Eu atoms. The maroon marked area corresponds 
to the surface-projected bulk band structure whereas the grey lines 
are the result of a calculation for a slab model. The green (blue) 
dots scale with ratio of localization: if the charge density is mainly 
located at the surface, then the size of the respective dot representing 
the eigenfunction is at its maximum whereas for equally distributed 
charge-density it vanishes. 



by the linearly dispersing valence band. Where "heavy" and 
"light" electron bands would cross, the energy degeneracy is 
lifted by hybridization leading to the formation of small en- 
ergy gaps and a mixing of the Eu 4f and Rh 4d states that 
is explicitly reflected in the ARPES data. The linear bands 
clearly hybridize with all J-terms of the Eu 4f multiplet. The 
resulting composite states exhibit a non- vanishing / charac- 
ter and are, thus, nicely seen at the chosen photon energy of 
120 eV where 4f emissions dominate. Effectively massless 
quasiparticles with finite 4f character appear in this way even 
at the Fermi energy and form a strong contrast to the heavy 
4/-derived quasiparticles usually observed in Kondo systems. 

Looking closely at the f-d hybrid band, running towards 
Ey in a quasi-conical shape, we can make important observa- 
tions. The energy position of the cone apex depends on the 
surface termination: for Si termination it appears « 60 meV 
above the Fermi level while for Eu termination it is located 
^65 meV below. Upon changing the photon energy no no- 
table dispersion in k z -direction is observed. These properties 
are well reproduced in band structure calculations. To this end 
we used a slab geometry where 55 (57) atomic layers parallel 
to the surface for Si (Eu) termination are considered. That 
allows to discriminate between surface and bulk effects by 
the localization of charge density. Treating the 4f electrons 
as core states the linearly dispersing bands are reproduced by 
electron states of almost two-dimensional character which are 
located at the surface. 

Understanding the interaction between localized 4f and 
itinerant electrons requires information on the strength of the 
hybridization which is intimately related to the symmetry of 
coupling orbitals. Within a simple hybridization model one 
may assume that the hybridization strength is proportional to 
the /-character of the unhybridized valence states at the site 
of the Eu atoms. 
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FIG. 3: LDA derived bulk band structure parallel to the re- 
direction for several planes k x x k y defined by k z . There are few 
different bands in the region around the T-point close to Ef. Going 
from the middle of the BZ to the top (bottom), one notices that they 
form a quasi-linear dispersion (2D) at k z = ±0.285 • 7r/a, which 
gets enhanced as well as shifted in energy at the surface. 



Respective information may be obtained from our slab cal- 
culations. It can be used then to simulate the ARPES spectra 
for both surface terminations of the EuRh 2 Si2 crystal. The 
results are shown in figs. [I]: and[T]l which qualitatively repro- 
duce the key features of our experimental ARPES data. The 
linear dispersion of the Rh bands obtained in the slab calcula- 
tions are mainly preserved in the hybridization model, modi- 
fied by hybridization gaps and the crossover to the Eu 4f 6 fi- 
nal state multiplet. Additionally, the hybridization model 
yields the transfer of Eu 4f spectral weight to the conically 
shaped valence band confirming the conclusions drawn from 
the experimental data. 

In fig. [2] we show the unhybridized bulk and surface bands 
for both surface terminations as computed from the slab 
model. In the calculations the localized 4/ states were treated 
as core states. Green and blue dotted lines represent surface 
bands, which are shifted in energy with respect to the bulk 
bands, and surface states, which are located inside a bulk band 
gap. The thicknesses of these lines scale with the surface to 
bulk ratio. The bulk band-structure projected into the surface 
BZ ifTTl is plotted in shades of maroon. There is remarkable 
agreement between the experimentally derived electron bands 
and the computed band dispersions for both terminations. In 
particular, the surface bands with linear dispersion at the T- 
point are well reproduced, including the energy difference of 
the cone apex position between the Eu and Si termination. In 
addition, there is a finite gap of 200 meV between the upper 
and the lower cone for both surface configurations. 

Group-theory considerations on two-dimensional lattices 
show that the quasiparticle dispersion around high- symmetry 
points is mainly determined by the symmetry of the space 
group of the crystal lattice lfl8l [T9l . Similar predictions for 
the space groups in three dimensions have not been carried 
out yet, but the arrangement of Rh atoms at the surface in the 
layered structure of EuRh 2 Si2 with tetragonal symmetry can 
up to some extent be approximated by a squared lattice for 
which in fact the appearance of linear electron dispersion was 
predicted. 

The layered character of the compound poses the question 
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FIG. 4: Experimental ARPES data of EuRli2Si2 recorded using 
950 eV photons, which allows to suppress the emission from the 
surface states and thus to disentangle the bulk electron structure of 
the system. The dotted rectangular shows the region where conically 
shaped Rh bands with their apex close to the Fermi level are coupled 
with the bunch of "horizontal" bands from the bulk Eu 4f 6 multiplet. 

whether similar phenomena as observed in the surface region 
may also be expected for the bulk. The projected bulk band 
structure shows indeed a band that coincides roughly with the 
linear surface band and reveals a weakly three-dimensional 
(3D) character (dark maroon color in fig. [2]). This becomes 
more evident in fig. [3j where we show cuts at several k z 
through the bulk BZ, parallel to the T-X-direction. For 
k z = 0.285 7i7a the bands in the vicinity of the r -point be- 
come degenerate and reveal linear dispersion while for other 
values of k z the band dispersions deviate slightly from this 
linear slope. Thus, similar hybridization phenomena as ob- 
served in the surface region may in fact be found in the bulk, 
too. To prove this, we performed an ARPES experiment in 
the soft X-ray photon energy range where the photoelectron 
escape depth is increased and the extreme surface sensitivity 
of XUV ARPES is partly overcome. 

In fig. [4] we show an ARPES map of EuRli2Si2 taken at 
his = 950 eV that extends over the T-points of three adja- 
cent Brillouin zones. As we know from respective results for 
the isostructural compound YbRli2Si2 l20l the contribution 
of the first subsurface Eu layer amounts to about only 20% 
of the total 4f emission at this photon energy range while it 
amounts to about 60% in the data shown in fig. [T] In fact, 
scanning the beam across the whole surface of freshly cleaved 
EuRfi2Si2 we did not detect a sign of surface-related electron 
states. In particular, we never found the surface Eu 4f feature 
at « 1.6 eV BE although the atomic cross section for Rh 4d 
and Eu 4f is of the same magnitude at 950 eV photon en- 
ergy. In the XUV photon energy range it was always seen for 
some regions of the sample surface. The absence of surface 
related features confirms the bulk sensitivity in the soft X-ray 
range. In the region around the T-points we observe well- 
defined bands that approach in a quasi-linear way similar 
to the surface band discussed above. Furthermore, we also 
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detect the string-like Eu 4f 6 final-state multiplet which again 
hybridizes with the linearly dispersive bands. The strong sim- 
ilarities between the surface- sensitive XUV and the more bulk 
sensitive soft X-ray data is owed to the layered structured of 
the material and the weak three-dimensional character of the 
underlying Rh 4d band. Although the energy and momentum 
resolution in the soft X-ray range are not as high as in the 
XUV they still allow the conclusion that ultra-light and heavy 
quasiparticles coexist for particular cuts through the 3D Bril- 
louin zone of the material, too, and are subject to the same 
hybridization phenomena as their surface counterparts. 

In summary, we have demonstrated that massless and heavy 
fermion quasiparticles may not only coexist but can also 
strongly interact in solids. Applying XUV and soft X-ray 
ARPES on EuRh 2 Si2 single crystals we observed linearly dis- 
persing Rh 4d-derived bands, corresponding to massless Dirac 
quasiparticles. Their conical shape is mainly imposed by the 
in-plane square symmetry of the layered compound. These 
linear bands are interfered by flat, dispersionless Eu 4f states 
corresponding to quasiparticles with extremely large effective 
masses. At the regions of intersection of both kinds of bands 
hybridization gaps as well as strong admixture of the 4f elec- 
tron states to the Dirac fermions are detected. This unambigu- 
ously demonstrates a strong interplay between these two op- 
posing limits of mass-renormalized fermions. Similar obser- 
vations as for the surface were made for the bulk of the crystal. 
Since the symmetry of the layered (122) compounds dictates 
the formation of ultra-light quasiparticles, one may anticipate 
that details of their interactions with 4f states could play a 
role in the emergence of unusual thermodynamic properties 



like Kondo physics or quantum criticality in other compounds 
of this family. 
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